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Abstract 



We present a self-contained theory for the mechanical response of DNA in single molecule 
experiments. Our model is based on a ID continuum description of the DNA molecule 
and accounts both for its elasticity and for DNA-DNA electrostatic interactions. We 
consider the classical loading geometry used in experiments where one end of the molecule 
is attached to a substrate and the other one is pulled by a tensile force and twisted by 
a given number of turns. We focus on configurations relevant to the limit of a large 
number of turns, which are made up of two phases, one with linear DNA and the other 
one with superhelical DNA. The model takes into account thermal fluctuations in the 
linear phase and electrostatic interactions in the superhelical phase. The values of the 
torsional stress, of the supercoiling radius and angle, and key features of the experimental 
extension-rotation curves, namely the slope of the linear region and thermal buckling 
threshold, are predicted. They are found in good agreement with experimental data. 
Key words: elasticity; self-contact; twist-storing polymer; DNA electrostatics 
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I Introduction 

Mechanics of the DNA molecule plays a key role in several biological processes at the 
cellular level. In several cases, the action of enzymes and proteins on DNA has been 
found to depend on the mechanical stress present in the molecule. For instance, the 
torsional moment in DNA controls the action of topoisomerases or RNA-polymerases 
[3). In this context, experiments where forces and torques are applied to a single DNA 
molecule provide a remarkable opportunity to gain insights into the mechanics of DNA. 
We are here interested in extension-rotation experiments using either optical or magnetic 
tweezers (0, 13) S 0, 0) j see Ref. (0) for a review. These experimental setups are equivalent 
from a mechanical perspective: a dsDNA molecule is fixed at one end on a glass pane 
while the other end is attached to a bead that pulls and twists it. In these experiments, 
traction and rotation are controlled differently: for the rotation mode, the twist angle 
is prescribed and the twist moment varies accordingly; for the stretching mode, the 
extension can vary although the pulling force is prescribed. DNA is under or over-wound 
and various molecule conformations are observed (0). In the present study, we focus on 
the over-winding of a dsDNA molecule under large imposed rotations: the molecule coils 
around itself in a helical way and forms plectonemes, as sketched in Fig. [H An important 
feature of the experimental loading curves is the linear decrease of the vertical extension 
of the molecule as a function of the imposed rotation. We have shown in previous 
studies [lol ) that this behavior can be captured by a purely elastic rod model based 
on Kirchhoff-Love elastic rod theory. In this paper, we extend this model and investigate 
the combined effects of elasticity and electrostatics. 

The response of plectonemic DNA under stress involves various physical phenom- 
ena such as elastic deformations, thermal fluctuations, electrostatic interactions, self- 
avoidance. Although some of these effects have been considered in the literature, a 
model addressing them all together is still lacking. Mechanical models of twisted rods in 



contact have been introduced, from an analytical (jlll ) or numerical ()12l ) perspective, but 
thermal fluctuations are not treated. A simplified analytical model, including some ac- 
count for fluctuations but omitting contact forces in the plectonemic region, is proposed 



in Ref. (Il3l ). Statistical mechanics of p lectonemic DNA has been approached, either 



analytically (|14| . Il5l ) or numerically using a Monte-Carlo method. The validity of 



some of these results was questioned in Ref. (|l7l . Il8l): in addition, long-range potentials 
raise convergence issues that have not yet been overcome in Monte-Carlo simulations. 
A composite model, gathering results from torsionally constrained polymer (flil ) and 
Monte-Carlo simulations has recently been introduced (jiol ). To date, this is the only 
model that confronts its predictions against experimental data. However it relies on an 



assumption on the supercoiling free energy that is not always valid ([211, l22l ) and uses 
parameters extracted from Monte-Carlo simulations. 

In this paper, we present a self-contained analytical model for the mechanical re- 
sponse of plectonemic DNA in extension-rotation experiments, which builds up on pre- 
vious work (j^). We focus on the plectonemic regime at large imposed rotations. This 
corresponds to the linear region in the experimental extension-rotation curves. Our elas- 
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Figure 1: Sketch of the experimental setup: a dsDNA molecule is fixed by one end 
to a glass pane while the other end is attached to a mechanical system, symbolized by 
the grey disc, which allows one to exert a pulling force F^xt and impose a rotation Inn. 
For large numbers of turns the molecule coils around itself in a helical way and forms 
plectonemes. The configuration of the molecule is made of two phases: the tails and 
the plectonemes. The plectonemic phase is characterized by superhelical radius R and 
angle a between the tangent t(s) and the helices axis. The dashed parts represent the 
regions we neglect: the matching region between the tails and the plectonemes and the 
end loop. 

tic model accounts for DNA-DNA interactions in the plectonemic region and for thermal 
fluctuations in the tail regions, where they are dominant. It captures the main features 
of the experimental curves and allows quantitative comparison to experiments with no 
adjustable parameter. 

The paper is organized as follows. In the next section we present our model and 
derive the equilibrium equations for a DNA molecule comprising plectonemes, for a 
generic interaction potential. In section IIIIl we describe two representative DNA-DNA 
interactions potentials available in the literature, which we then plug into our model. 
The results are then compared with experimental data in section IIVI 
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II Model 



Our description of the DNA molecule is based on a coarse-grained representation (? ). 
We introduce a continuum rod model whose mechanical behavior is similar to that of 
the molecule, and makes use of effective elastic and electrostatic properties obtained by 
smoothing out the details at a scale of several base pairs. We deal with an inextensible 
elastic rod with circular cross section, bending rigidity Kq, and twisting rigidity K-^. 
The loading geometry is that of Fig. [H and applies to the experiments where the lower 
end of the molecule is clamped on a glass pane and the other end is subjected to a 
tensile force Fext and rotated by n turns {i.e. an angle 27rn). The imposed rotation is 
achieved through a torsional moment Mgxt- Note that the torsional moment has become 
accessible to experimental measurements (0) only recently. 

Geometry 

The inextensible rod, of length i, is parametrized by its arc-length s, the origin s = 
being at the lower end. The rod centerline is described by a vector-valued function 
r{s) and its unit tangent t{s) '= dr/ds. The geometric curvature of the rod is noted 

k(s) |dt/ds|. The twist is noted t(s): it describes the relative rotation of neighboring 
cross sections about the tangent t{s). Note that the twist is a different quantity from 
the Frenet (geometric) torsion of space curves — the latter is irrelevant in the context 
of elastic rods. 

We consider the geometry of the double stranded DNA sketched in Fig. [H which 
is relevant to the plectonemic regime: two twisted, straight tails are separated by a 
plectonemic region composed of two identical and uniform helices. Note that each helix 
is itself a piece of the double stranded DNA molecule. For a large number of turns n, 
the loop at the end of the plectonemes and the curved region connecting the tails to 
the plectonemes are much smaller than the tails and the helical parts, and hence are 
neglected. Even though we depict the plectonemic region as a single chunk for simplicity, 
our model applies equally well to the case where the plectonemes are distributed in 
several places along the molecule; then, the elastic rod is made up of two phases, one 
with linear DNA and the other one with plectonemes, and the plectonemic structure in 
Fig. [T] represents the plectonemic phase collectively. The molecule contour length spent 
in the tails phase and in the plectonemes are noted it and ip respectively. They sum up 
to the total length i = ip + it. The plectonemic phase is characterized by its superhelical 
radius R and its superhelical angle a, which are assumed to be uniform: neither R nor a 
may depend on s, although they depend on the loading. Curvature is zero in the straight 
tails, and takes a constant value in the plectonemes which can be evaluated using simple 
geometry. The integrated squared curvature, which enters into the bending energy, is 
then found to be: 



see Ref. ((91). Since the rod has a circular cross section, the twist t{s) is uniform, dr/ds = 




(1) 
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0, as shown for instance in Ref. (|23l ). The internal torsional moment M{s) in the rod is 
related to the twist r(s) by the constitutive law M(s) = i^3r(s). Therefore, its value 
M(s) is constant along the rod, and equal to the torque M^xt = r applied by the 
bead. In what follows, we study the equilibrium of the rod and compute the parameters 
R, a, £p and r = Mext/-f^3 as a function of the loading (pulling force Fcxt and number 
of turns n) by minimizing the energy. 



Variational formulation 

We derive the energy of the system as a function of the superhelical angle a and radius 
R, of the twist r and of the plectonemic contour length ip. Equilibrium solutions and 
their stability will be derived later on by minimizing this energy. The experiments are 
performed under imposed end rotation: energy minimization is performed under the 
constraint that the number of turns n imposed on the bead is equal to the link Lk of the 
DNA molecule. Neglecting the writhe of the tails, the link can be written as 0,0): 

rr. 1 /"^ 1 slu 2a 1 / „ slu 2q „ \ 

„ = Lk = T„ + Wr = -y__.d.-x^ = -(rf-x^f„j • 

where x = il stands for the chirality of the two helices of the plectonemic phase. 

The total energy of the system is the sum of three terms, V = Vei + Vext + ^mtj where 
the first is the strain elastic energy, the second is the potential energy associated with 
the external load -Fextj and the last term accounts for DNA-DNA interactions between 
the two helices in the plectonemic phase. The strain elastic energy of the rod is defined 
as the sum of a bending term, proportional to the integrated squared curvature, and a 
stretching term, proportional to the integrated squared twist: 

T/ ^0 f'^ 2. , 2. i^'osin'^a. , ^ 

This energy captures the elastic behavior of the rod in response to applied forces and 
moments; it is zero in the natural (straight, twistless) configuration of the rod. The 
pulling force is described using a potential energy: 

^cxt = -i^cxt {Z{£) - Z{0)) = -Fext Az . (4) 

Here Az '= {z{i) — z{0)) is the extension of the molecule along the direction z of 
application of the pulling force. Since we assume the tails to be straight and neglect the 
curved region region connecting the tails and the plectonemes, the vertical extension of 
the filament reads Az = it = i — ip and we can rewrite: 

Vext = -Fext {i - ip) . (5) 

There is no need to consider any potential energy associated with the rotation of the 
end attached to the bead since the energy will be minimized for a given rotation of the 
bead using the constraint on the link. 
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In previous work ((91) we solved this elastic rod model by assuming the superhelical 
radius to be prescribed and extracted its value from experimental data. Here, we take 
a more principled approach and complement the above elastic equations with a proper 
model for DNA-DNA interactions in the plectonemes; in particular, this makes it possible 
to predict the superhelical radius. These interactions are dominated by different physical 
effects depending on the separation distance between the two DNA superhelices. In 
the range of separations relevant to extension-rotation experiments, of order several 
nanometers, electrostatic effects dominate. In our model, interactions are limited to the 
plectonemic phase and are described by an energy contribution of the form: 

Mnt = ip U{R, a) . (6) 

This energy depends on the superhelical parameters R and a, and is proportional to the 
plectonemic contour length £p, and assumption valid when ip is much larger than R. 

The total energy of the system is defined as the sum of the elastic, potential and 
interaction contributions: 

V{R, a, £p,T) = ^^^p + ^r'^ - F,^t {i - + U{R, a) . (7) 

It will be minimized subjected to the end rotation constraint given by Eq. [2l This 
constraint provides an affine relation between n and ip and so can be used to eliminate 
the quantity ip in favor of n. Dropping the constant term (— -Fext i) in the energy, we 
obtain: 

V{a,R,T) = ^T^i + {2Trn-Ti) 



—2x f Kq sin a 
sin 2a I ^ R~ 



+ RFe,t + RU{R, a) 



(8) 



Equilibrium equations 

The total energy of the system, given by Eq. [HI takes into account the fixed end rotation 
since Eq. [2] has been used to eliminate ip. The equilibria of the rod are then directly 
given by minimization of V{a, R, r) with respect to its three arguments. Canceling the 
first variation of V, that is writing §^ = 0, = and ^ = 0, we obtain 

cosasin^Q dU(R,a) 2 / i^'n sin^ a „ „ x\ / ^ 

Fext - ^ sin^ a + R ^^^^'"^ + U{R, a) = , (9b) 

Mext + f + ^i^cxt + RU{R, a)]=0. (9c) 

sm 2a \ 2 R J 

In the first term of the last equation, we have eliminated r in favor of Mext using the 
constitutive relation Mext = r, and thereby removed the twist rigidity from the 
equations (its value is not known with good accuracy). 

The set of three nonlinear equations ([9]) must be solved for the three unknown values 
of the parameters q, R and Mext at equilibrium, given the value of the external force 
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Fext- This requires an interaction potential U{R, a) to be specified, as is done in the 
forthcoming Section IIIII The set of equations ([9]) extends the model of Ref . (0) , vahd 
for non-penetrable tubes, to filaments in long-range interaction (such as electrostatic 
interactions) . 

Note that the equations ([9]) do not depend on the number of turns n. As a result, 
their solution a, R, Mgxt do not depend on n either. The equations describes the 
equilibrium of two phases; increasing n transfers some arc length from the tail phase to 
the plectonemic phase, without changing their properties. This invariance with respect 
to n explains the presence of a linear region in the experimental curves, as shown in the 
next section. 

The term in parenthesis in Eq. [9c] is always positive. This shows that the sign of 
the chirality x = =tl is opposite to that of n: rotating the bead in the positive direction 
n > for instance, requires a positive torque Mext, hence a negative x = — 1 by this 
equation (left-handed super helices). 



Vertical extension of the filament 

In extension-rotation experiments the vertical extension of the filament is recorded while 
the number of turns is increased. The formula Az = it, valid for straight tails, does 
not holds in the presence of thermal fluctuations. We account for these fluctuations by 
introducing a rescaled quantity: 

Azth = Pwic Az , (10) 



where the correcting factor p^ic is given by the worm-like chain theory ([2^) as the solution 
of: 

( A:bT)2 - ^-^^ + 4 (1 - 4- ^''> 

Here A:^ is the Boltzmann constant and T the absolute temperature. In order to write 
Azth as a function of the number of turns n, we use the equality Az = it = i — ip in the 
right-hand side of Eq. \lU[ and use for ip the expression obtained by solving Eq. [2j 

/ 2R \ ^ 4ttR , , 

Azth = 1 - X ■ o /^wlc i^X Pwlc ■ o ■ (12) 

\ sm 2a / sm 2a 

Recall that neither a, R nor r = Mext/-f^3 depend on n; as a result, the extension Azth 
depends linearly on the number of turns n in the above equation. This linear dependence 
is a well-known feature of the experimental curves. 



Ill DNA-DNA interactions 

In the variational formulation exposed in the previous section we have introduced an 
energy U{R, a) describing DNA-DNA interactions. At moderate distances DNA-DNA 
interactions in solution mainly originate from electrostatic effects between the charged 
sites of the two molecules (phosphate groups) and between these charged sites and the 
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counter and co-ions present in the solution. The theoretical analysis of the long and short 
range interactions between two poly-ions in solution has been the subject of numerous 
studies (2^, H) and there is currently a regain of interest in this question due to recent 



progress in single-molecule experiments — see (j27l ) for a review. In the present model 
the interaction energy U{R, a) is specified independently of the mechanical behavior of 
the molecule. As a result, we can combine the elastic description of the previous sections 
with different theories for DNA-DNA interactions. In the following, we illustrate this 
approach using two representative interaction energies U{a,R) that can be found in the 
literature. 

We favor interaction energies U{R, a) that can be expressed in closed analytical form 
and have no adjustable parameters; this enables us to make predictions and compare 
them to experiments, rather than to fit existing data. In the literature on DNA-DNA 
interactions (27, 2^, 2^), we picked two well-established models satisfying those require- 



ments. The first one, UpB{R,Oi), derives from the Poisson-Boltzmann equation and 
was obtained by Ubbink and Odijk (jiol ): the second one, Ucc{R^Oi), is based on the 



counterion condensation theory (j3ll ) and was derived by Ray and Manning (j32l ). These 
two models address the electrostatics of DNA in solution but their treatment of the 
interactions between DNA and the ions in solution differ substantially. 

Poisson-Boltzmann model 

In their study of supercoiled DNA plasmids (j3Q ) . Ubbink and Odijk derive an analyt- 
ical expression for the electrostatic interaction energy between two interwound DNA 
molecules. Their work is based on the Poisson-Boltzmann framework (PB); in the com- 
putation of the electrostatic repulsion of the two charged molecules, the presence of the 
counter-ions and co-ions in solution is considered. It has been shown in Ref. (jsil ) that 
the non-linear PB problem could be simplified to a linear one by considering screened 
(Debye-Huckel like) potentials and renormalized linear charge densities v. The value of 
the effective charge v is obtained by matching the solution of the non-linear PB equation 
with the solution of the linear PB equation in the far-field region. However there is no 



consensus on the exact value of this effective charge and the various numerical ([SJ, l35l ) 
or analytic al ( 361 . 1371) studies yield scattered results. 

In Ref. (1301 ) ■ Ubbink and Odijk compute the electrostatic interaction energy per unit 
length as 

UpB{R,a) = hBTu^lB^^^e-^^^^'^^^{a) , (13a) 

where the angular dependence reads 

ip{a) = 1 + 0.83 tan^ a + 0.86 tan^ a . (13b) 

Here, kB is the Boltzmann constant, T the temperature in Kelvin, the effective linear 
charge density (in m~^), Ib the Bjerrum length, and the Debye length. For a typical 
temperature T = 300 K we have Ib = 0.7 nm, and for a monovalent salt concentration 
c = 10 mM the Debye length is = 3.07 nm. The value of the effective charge f 



Elasticity and electrostatics of plectonemic DNA 



9 



depends on salt concentration, its value for a monovalent salt concentration c = 10 mM 
is taken as = 1.97 nm~^ according to Ref. (jsol ). 

The calculation of the interaction energy can be simplified by taking a = 0, hence 
(p{a) = ip{0) = 1, which amounts to consider two straight and parallel molecules; this 
approximation has been used for instance in Ref. (0). In the rest of the paper Up^ will 
refer to the potential obtained under this approximation, namely U'p^{R) = Upb{R,0)- 

Ray and Manning model 

The interaction energy derived by Ray and Manning in (js^ ) is based on the counter-ions 
condensation theory (|3ll). The authors examine the interaction of two straight and par- 
allel DNA molecules with charged sites in solution (the dependence on the superhelical 
angle a is neglected). The main point of the theory is to consider that part of the DNA 
bare charge is neutralized by the condensation of the counter-ions around the molecule. 
The energy is the sum of three terms: interactions between pairs of charged sites belong- 
ing to the same DNA segment, interactions between pairs of charged sides located on 
opposite segments, and adsorption energy of the counter-ions onto the molecule. Three 
cases are considered, namely short, intermediate, and long interaxial distances between 
the molecules. The short distance case, below the crystallographic radius of DNA, is 
not relevant to our analysis. The intermediate case introduces an adjustable parame- 
ter, which we try to avoid. Consequently we only use the long-range case, relevant for 
inter-distances larger than the Debye length; in our notations it writes: 



where 6 = 0.17 nm is the charge spacing parameter of the DNA molecule, and i = Is/b 
is the dimensionless charge density of DNA (.^ = 4.11 at T = 300 K). The function 
Bj^{x) is the modified Bessel function of the second kind and order 0. 



We solve equations [9] for the superhelical radius R, angle a, and external torque Mext, 
using one of the interaction energies Upb, Up^ or Ucc- These equations are nonlinear 
and their roots are found numerically using a Newton-Raphson algorithm. We present 
the results for the superhelical variables R and a in Fig. [21 for the torsional moment 
Mext in Fig- El We also plot derived quantities, to be defined later, such as the slope 
q of the extension-rotation curves in Fig. IH and the thermal buckling threshold n* in 
Fig. El We compare our results with the model of Ref. (fiol ) and with experimental data. 
To ease comparison with our previous work (0), we use the same set of experimental 
data. These data were obtained on a 48 kbp lambda phage DNA molecule in a 10 mM 
phosphate buffer. 

With the interaction energies used in the present paper, we find that the nonlinear 
equations have two roots below a threshold value of the force, and no root above. For 




(14) 



IV Results 
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Figure 2: Computed values of the superhelical radius R and angle a (inset) as functions 
of the pulhng force Fext, using one of the interaction energies Upb, Up^ or Ucc- These 
plots are obtained by solving the equilibrium equations [9] for each value of the pulling 
force -Fext- The function a(-Fcxt) decreases for all interaction energies at large enough 
forces — with the energy Ucc^ this decrease occurs beyond the domain of forces shown 
in the figure. 



a salt concentration c = 10 niM, this threshold value of the force is 4.7 pN using Upb, 
4.9 pN using Up^, and 6.9 pN using Ucc', all these values are above the maximum pulling 
force applied in typical experiments. We have studied the stability of the two solutions 
corresponding to the two roots of our equations when the force is below threshold, and 
found that that with lower a and R is unstable; the other one is stable. We study and 
plot the latter in the following. 

When the force reaches its threshold value, an instability occurs and the stable 
solution disappears by merging with the unstable one. For larger forces, no stable 
solution exist and the two helical parts collapse. This may be related to the observation of 
tightly supercoiled configurations in Ref. (jSa ). The collapse arises when the electrostatic 
interaction is no longer strong enough to sustain the applied force; the description of 
collapsed solutions would require an account of hard-core repulsion and other short-range 
forces. 

Superhelical geometry 

The quantities R and a are plotted Fig. [2] as a function of the applied force Fcxt • The 
curves R{Fext) and a(-Fcxt) obtained for the different interaction energies Upb, Upp or 
Ucc are close over the entire range of forces. As will be confirmed later, the predictions 
based on the different interaction models are all very similar. 
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Figure 3: Computed values of the torsional moment in the molecule Mcxt as a function of 
the pulling force Fcxt • We compare the results of the present model using the interaction 
potentials Upb, Up^ or Ucc, to predictions in Ref. (0), namely Marko's Eq. 17 with 
^ = 50 nm, C = 95 nm and P = 28 nm, and Strick et al. (|4ll) . and results from our 
previous work (0). 



As expected, the superhelical radius decreases with the pulling force; note that it 
becomes less than the Debye length for forces above F^xt ~ 2.5 pN. The superhelical 



angle a is known to be a control parameter in the action of the topoisomerases (l39l'l. I t 
is plotted in the inset of Fig. [2j In contrast with models of elastic tubes in contact (|4d). 
where a increases monotonically and reaches the value 7r/4 asymptotically at large forces, 
we find here that it reaches a maximum well below 7r/4 and then decreases, due to long- 
range forces. This decrease has already been observed in Ref. (0), where the value of 
the superhelical angle was extracted from experimental data. 

Torsional moment 

Recall that the torque Mext = r applied by the bead in order to impose a rotation 
27r n does not depend on the number of turns n by equations. [9l This quantity is plotted 
in Fig. [3]as a function of the pulling force -Fcxt- We compare our results (i) with Marko's 



Eq. 17 in Ref. (j2Cll ): (ii) with a formula Mext = \/ 2 Kp F[Zf, obtained by Strick et al. by 



approximating the plectonemes by a chain of circles (l4ll ) : and (iii) with ourprevious 
study (jH) based on hard- wall interactions. Even though Refs. (|2ol). (41) and Q do not 
address long-range interactions, all the curves reveal a similar behavior: the moment 
increases monotonically with the applied force, with a decreasing slope. However, our 
results show that long-range interactions significantly increase the value of the moment 
required to achieve a given rotation. 
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Figure 4: Computed values of the slope of the plectonemic region g as a function of 
the pulling force -Foxt- Experimental points with filled circles are extracted from the 
experimental hat curves shown in Fig. 2 of ([^). The thin purple curve is obtained from 
the theory in (fiol ). using the same parameter values as in Fig. [3l 



Extension-rotation curve 

Our model predicts that the derivative of the vertical extension Azth with respect to 
n is constant, i. e. that the extension-rotation curve is linear in the regime of large 
rotations that we consider. This linear regime is well-known experimentally, see Fig. [SJ 
From Eq. [121 t^ie slope q is given by: 

dAzth 47ri? 

q = —5 — = x^^/5wic ■ 15 

an sm 2a 

Its value is computed using the values of a and R obtained earlier by solving the equi- 
librium equations. We plot in Fig. H] the slope q as a function of the force, for the 
three interaction energies. For comparison, we also plot the slope predicted by Marko's 
model (0), and the slopes read off directly from experimental extension-rotation curves 
(these experimental data were kindly provided by V. Croquette and have appeared in 
Fig. 2 of our Ref. 0) and in Refs. (0) and 0)). 

Our model shows good agreement with the experimental data which are reproduced 
in a more consistent manner, especially at low forces, than in Ref. (0). In this reference, 
hard-wall interactions are used with a supercoiling radius independent of the pulling 
force; this may be the cause of the poorer agreement with experimental data at low 
forces, when long-range interactions dominate. 

A typical extension-rotation curve comprises two regions: a linear region for large 
n, which we have been discussing so far, and a parabolic region at low n studied in 
Ref. (ji^). The central region is dominated by thermal effects and will not be addressed 
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Figure 5: Experimental curve showing the vertical extension Az^xp of a lambda phage 
DNA 48kbp molecule as a function of the imposed number of turns n, at constant force 
-^ext = 0.44pN. The quantity Azth{n) defined in Eq. [T2]is our prediction for the linear 
part of the experimental curve. The number of turns at the transition n* and the slope 
q are also shown. 



here. However, we can characterize the transition between the two regions. The number 
n* of turns at which the transition occurs is defined using the linear extrapolations shown 
in Fig. El This n* corresponds to the onset of the plectonemic regime. In our model, 
it is computed by setting £p = in Eq. [2l this yields n* = t£/{2tt). Recall that the 
value of r = Mg^t / is computed from the equilibrium equations [9l We plot in Fig. [6] 
the value of n* as a function of the force and compare to the values extracted from the 
experimental curve as well as the value from the theory in Ref. (0). 



V Discussion and conclusion 



In Section IIVI we found that the solutions disappear above a threshold value of the 
pulling force. This can be interpreted as the fact that the physical solution involves 
a very short inter-distance, although we have retained the long-distance part of Ray 
and Manning's potential only. This can be cured in principle by restoring a complete 
expression of the potential given in Ref. (32). In its complete form, the potential is 



non-monotonous and several inter-distances are possible for a given value of the control 
parameter. This feature indicates the possibility of a transition from a classical super- 
coiled state to a tight super coiled state. A possibly related transition has been reported 



^Except in the long-distance part used here, Ray and Manning's potential involves physical quantities 
that are unknown. 
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Figure 6: Computed values of the number of turns n* at the transition as a function of 
the puhing force -Fextj using = 95 nm ksT. Experimental points are extracted from 
the curves shown in Fig. 2 of (0). The thin purple curve is obtained from the theory 
in liol ) , using the same parameter values as in Fig. O 



in experiments (jSa ). When used in conjunction with the complete Ray and Manning's 
potential, our model could provide a bridge between the analytical expression for the 
potential and the experiments, and provide a quantitative account of the transition. 

We have presented an analytical model for DNA supercoiling in extension-rotation 
experiments. It is based on an elastic description of DNA deformations, carefully ac- 
counts for DNA-DNA interactions in the plectonemic region, and makes use of a valid 
formula for the link. DNA interactions are modelled using long-range forces computed 
from potentials available from the literature. Their description is compartmentalized 
from the rest of the theory, which makes it possible to test different interaction ener- 
gies. We have used our model in combination with two interaction energies. These 
energies come from different physical contexts and are widely used in the literature. 
Using either one, we find good agreement with experimental data without adjusting any 
parameter. This suggests (i) that using yet another energy for electrostatic interactions 
would yield comparable results and (ii) that DNA tertiary structures are determined to 
a large extent by the elasticity of the molecule and do not depend heavily on the details 
of the interaction model. The extension-rotation experiments that are now routinely 
performed can then be viewed as a way to probe the elastic properties of the molecule. 
Given that the mechanics of DNA under combined twist and tension can be captured 
by a relatively simple analytical model, an interesting direction for future research is to 
extend this model to the mechanical action of proteins, such as RecA for example, on 
super coiled DNA. 
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